Theory of extraordinary optical transmission through subwavelength hole arrays 



L. Martm-Morenoi, F. J. Garci'a-Vidal 2, H. J. Lezec^, K.M. Pellerin T. Thio ^, J. B. Pcndry 5, T. W. Ebbesen ^ 
^ Departamento de Fisica de la Materia Condensada, ICMA-CSIC, Universidad de Zaragoza, E-50015 Zaragoza, Spain. 
^ Departamento de Fisica Teorica de la Materia Condensada, Universidad Autonoma de Madrid, E-28049 Madrid, Spain 

^ ISIS, Universite Louis Pasteur, 67000 Strasbourg, France 
NEC Research Institute, Princeton, New Yersey 08540, USA 
^ The Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 

We present a fully three-dimensional theoretical study of the extraordinary transmission of light 
through subwavelength hole arrays in optically thick metal films. Good agreement is obtained with 
experimental data. An analytical minimal model is also developed, which conclusively shows that 
the enhancement of transmission is due to tunneling through surface plasmons formed on each metal- 
dielectric interfaces. Different regimes of tunneling (resonant through a "surface plasmon molecule", 
or sequential through two isolated surface plasmons) are found depending on the geometrical pa- 
rameters defining the system. 

PACS numbers: 78.66.Bz, 73.20.Mf, 42.79.Dj, 7f.36.-)-c 



In the last few years, and mainly due to advances in 
nanotechnology, there has been a renewed interest in ex- 
ploiting the dielectric response of metals to make pho- 
tonic materials For instance, the photonic in- 
sulating properties of metals can be used to trap in- 
cident radiation, focusing light in very small volumes 
iQ-^ . Very recently , another interesting effect of light 
interacting with structured metals has been discovered: 
the transmission of light through subwavelength hole ar- 
rays made in a metal film can be orders of magnitude 
larger than expected from standard aperture theory 
Apart from its fundamental interest, this extraordinary 
transmission effect has potential applications in 
subwavelength photolithography, near-field microscopy, 
wavelength-tunable filters, optical modulators, and fiat- 
panel displays, amongst others. While the wavelength 
at which some transmission features appeared suggested 
that surface plasmons (SP) were involved in the 
process, the physical mechanism for the huge enhance- 
ment has not yet been elucidated. Some calculations 
have been performed for a simpler geometry: an array 
of slits [^ - |l5| , where high transmission was also pre- 
dicted. However, although interesting in their own right, 
these results do not apply to the experimental situation. 

In this letter, we present the first fully three- 
dimensional theoretical study of transmission through 
hole arrays, obtaining an excellent agreement with exper- 
imental data. Moreover, we develop a simplified version 
of the model that clearly captures the physics involved. 

Fig. 1 shows the experimental "zero-order" transmit- 
tance of light at normal incidence (Too), through an array 
of holes in a free-standing metal film. The free-standing 
metal film, of which the fabrication is described elsewhere 
consisted of a 220nm thick Ni core, perforated with a 
square array of holes by focused-ion beam milling. The 
film was subsequently overcoated with 50nm of Ag on 
both sides by sputter deposition which resulted in a coat- 
ing of the walls of the holes as well as the in-plane sur- 
faces of the film. The total thickness of the film was 



h — 320nm and the lattice constant of the hole array 
was L — 750nm. After coating the holes had an av- 
erage diameter of 280nm. It has been shown ||] that 
such a "sandwich" structure has the same transmission 
properties as an equivalent perforated film made of silver 
throughout. The advantage of a free-standing metal film 
is that it is possible to mill much better defined holes than 
in a film on a substrate. Moreover, in a free-standing film 
the dielectric constant is the same in all non-metallic re- 
gions. In Fig. 1, Too shows the well known Rayleigh min- 
ima of Wood's anomaly | ]l6t which appear in any diffrac- 
tive array roughly when an order of diffraction emerges 
tangent to the array. The data also shows the extraor- 
dinary transmittance effect: at A « SOOrtm, Too of the 
order of 15% was found. As the area covered by holes is 
only 11%, the normalised-to-area transmittance of light 
is 130%. Standard aperture theory for a single hole pre- 
dicts transmission efficiencies of order 1%. 

In our theoretical model, we consider a metal film of 
thickness h, perforated with square holes of side length d, 
periodically distributed, as in the experiment, in a square 
array of lattice parameter L. In the long- wavelength 
limit, experimental results have shown that the trans- 
mission coefficient depends on hole area, but does not 
appreciably depend on hole shape. The choice of square 
holes is, therefore, inessential, and was motivated purely 
for analytical convenience. When comparing with exper- 
imental results we take the square side as the square root 
of experimental hole area. We consider the general case in 
which light impinges from a uniform medium 1, into the 
perforated metal film (region 2) which lies on a substrate 
(medium 3), where the light in collected. The dielectric 
constants are, ei in medium 1, €2 inside the holes, and 63 
in medium 3. The metal is characterised by a frequency- 
dependent complex dielectric constant e„i(a;). 

We treat the EM fields in the metal through the surface 
impedance boundary condition |17|, a procedure fully 
justified, as the frequencies considered are well below the 
metal plasma frequency. We expand the EM fields in 
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the eigenmodes of the different regions. That is, the EM 
fields are a hnear expansion of S and P plane waves in 
regions 1 and 3, and Bloch waves combining (evanescent 
and, if there is any, propagating) TE and TM waveguide 
modes inside the holes. Although, in principle, an 
infinite number of modes in each region should be taken 
into account, results quickly converge with the number 
of parallel wavevector components considered. The prob- 
ability amplitude, Uf, for an arbitrary incident (from 
medium f) EM plane wave, i, to be transmitted to a 
outgoing plane wave /, in medium 3, was calculated in a 
multiple scattering formalism. By appropriately match- 
ing the EM fields in the 1-2 and 2-3 interfaces, tif can be 
expressed in terms of transmission and reflection ampli- 
tudes for a single interface as: 
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Latin indexes refer to plane waves, either in medium 1 
or 3, and Greek indexes to waveguide modes inside the 
holes. T^'^, T^^,p^ and are components of the scatter- 
ing matrix for a single interface, either 1-2 or 2-3. More 
precisely, r^^ jg g, matrix giving the transmission ampli- 
tude between modes in medium 1 to modes in medium 2. 
p^^ is the reflection amplitude for the waveguide mode 
a, traveling towards medium 1, to be reflected into the 
mode 7 traveling away from medium 1, after scattering 



with the 2-1 interface. 
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the interface 2-3. rjy connects mode (3 with outgoing 
mode /, through interface 2-3. Finally, = 6*''=°'', z is 
the direction perpendicular to the interfaces, and qza is 
the wavevector component along the z-direction for the 
waveguide mode a. 

Figure 2 shows the computed results for Too- Calcu- 
lations were done for the nominal experimental values 
for the different parameters defining the structure. The 
metal dielectric constant has been taken as the measured 
one for Ag |18|. It must be stressed that in our calcu- 
lation there are no fitting parameters. Also shown are 
results for slightly different values for hole size, d, in or- 
der to give an idea of the sensitivity of the results to 
non-uniformities of hole area. As shown in Fig. 2, our 
calculation also shows peaks, with Too which exceeds by 
2-3 orders of magnitude the predictions by a theory based 
on independent holes |9). Moreover, the wavelenghts at 
which peaks occur closely match the experimental values. 
Clearly the mechanism responsible for the extraordinary 
transmission is present in our model. Discrepancies be- 
tween calculated and experimental values for linewidths 
and maximum transmittances, can well be due to small 
variations in the hole diameter throughout the sample, 
as supported by the sensitivity of the transmittance to d 
shown in Fig. 2. 

Although the calculation also shows extraordinary 
transmission through periodic subwavelength holes in a 



metal, it is not obvious what mechanism is involved. 
However, as we show in the following, a highly simplified 
version of the model, which can be analytically worked 
out, also shows the same behaviour in the regime of in- 
terest: A > L S> d. In this minimal model, a strong 
truncation in parallel wavevectors is made: in regions 
1 and 3, only first order diffraction is considered, i.e., 
the possible wavevectors in the direction of the incoming 
electric field, x, are k^o — and kx±i — ±X' region 
2, where all modes are evanescent for L > 2d, only the 
most slowly decaying evanescent mode, the TEqi mode, 
is taken into account. Figure 3 shows the comparison 
of the fully converged calculation (solid line) with the 
results obtained with the minimal model (dashed line), 
for the nominal parameters of the experimental set up, 
both for (3a) the experimental dielectric constant for Ag 
and (3b) for an hypothetical case of non-absorbing Ag 
(/m [e((jj)] = 0). The minimal model captures well the 
extraordinary transmission phenomena: the main effect 
of truncating the number of modes is a small unimpor- 
tant shift in the position of the transmission peaks. In or- 
der to further simplify the discussion, in the following we 
consider the ideal non-absorbing case. Finite absorption 
merely reduces the height of the peaks (for zero absorp- 
tion the transmittance can be as large as 100 %) without 
altering the physical picture. 

For the symmetric case, ei = £2 = £3 = f , Pn — Pii = 
p and too simplifies to 
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where r^f = 2^o/(G2 + Gi), rfg = 2G2/(G2 + Gi) and 
p = (G2 — Gi)(G2 + Gi). G2 and Gi are "conductances" 
(inverse of total impedances) in region 2 and in region 
1, respectively. They have somewhat involved analytical 
expressions that, for the case of a hole array in a perfect 
metal, simpHfy to G2 — Qzi/g, Gi = Sq + 2Slg/kzi, 
with g = uj/c,So = d/L,Si = Sosiii[kxid/2] /{kxid/2), 
and kzi = 

Wood's anomaly is present in this model, associated to 
a zero in either rgf or r^g , i.e., when the conductance in 
medium 1 or 3, respectively, is infinite. 

With respect to the extraordinary transmission, trans- 
mittance peaks « 1 seem to point to the presence of 
resonant phenomena. This idea is reinforced by the pres- 
ence of multiple scattering denominators in Eq.(2) that 
could be close to zero. Naively, one would expect that 
this possibility never occurs: e~2i<?2il^ ^ 1 in the regime 
considered and p is the amplitude for reflection, which 
can be expected to be always less than 1. Actually, this 
is not true, as the condition \p\ < 1 applies only for prop- 
agating modes, where it corresponds to current conserva- 
tion. In the regime we are considering, all modes inside 
the hole are evanescent and current conservation only re- 
stricts Im [{p] > 0, with no restrictions on the real part 



2 



or the modulus of p . This opens up the possibiHty of 
resonant denominators for |p| ^ 1 . In fact, it can be 
shown analytically that \p\ has a peak at A slightly larger 
than L, with a maximum value that scales as Fig- 
ure 4a shows IpI as a function of A for the parameters 
considered in Fig. 2. It is illuminating at this point to 
consider the dependence of the transmittance on h. Fig- 
ure 4b illustrates graphically that the peaks in Too occur 
at the A for which the distance between \p\ and e'''^^''' is 
minimal. For zero absorption, a Too close to unity occurs 
whenever this distance is zero, i.e. when the curves in- 
tersect. For large enough h, the curves do not intersect 
and Too maxima decay exponentially with h. 

All this is mathematics, showing that the effect of ex- 
traordinary transmission has a resonant nature. But, 
what is the physical origin? Analytically, we find that the 
frequency at which the maximum of |p| appears coincides 
with the SP frequency of a periodically (air-metal with 
holes) isolated interface at parallel momentum ^ (what 
we are going to call first-order SP). This is a localised 
EM mode, confined to the interface. It is not the only 
EM mode at that frequency, as there also exists the zero- 
parallel-momentum propagating mode, into which the SP 
can decay. This decaying channel provides a linewidth 
and a maximum value to p(A), which otherwise would be 
a delta function. Furthermore, when we calculate the SP 
of the periodically perforated metal slab, we find that the 
single interface first-order SPs on the two metal surfaces 
combine to form a " SP molecule" , in much the same way 
as electronic states of isolated atoms combine to form 
molecular levels. The frequencies of the first-order SP 
molecule are given precisely by the same condition for 
resonant denominators in the expression for too, if the 
coupling to the radiative modes is neglected. As in the 
case of the isolated interface, this SP molecule coexists 
with a continuum and, rather than an exact eigenstate, 
is a (more or less long lived) resonance. 

This point on the formation of the SP molecule can be 
more precisely discussed by considering typical times in 
the transmission process. The width of the peak of p{ijj) 
is related to the lifetime of the SP, t^ad, the time taken 
for a SP of an isolated surface to decay into the radiative 
mode. We obtain that, approximately, trad ~ (c|fc2i|)~^ 
. From the frequency splitting between molecule lev- 
els, taken as uncoupled from the radiation modes, the 
time needed for each of the molecular modes to form, 
tres, can be estimated as Ues ~ trad{d?'e^'''^^'^ / L'^)- 
So, in the small h regime, where the intersection of \p\ 
with the exponential occurs at the bottom of the \p\ peak, 
ires "C trad and the molecule levels are fully formed. The 
photon then goes back and forth several times inside the 
hole, building up coherent constructive interference in the 
forward direction much as would occur in electron reso- 
nant tunneling or in a Fabry-Perot interferometer (but 
with evanescent waves instead of propagating ones). As 
we consider larger h, tres increases, while trad, being a 



property of an isolated interface, essentially remains con- 
stant. So, the photon can make fewer round trips inside 
the hole before being radiated to infinity, and the con- 
cept of plasmon molecule becomes less well-defined. The 
condition tres — trad marks the value of h at which the 
transmittance maximum ceases to be one. For even larger 
h, tres ^ trad and the process is more like sequential tun- 
neling, where the incoming photon gets trapped in a SP, 
tunnels to the SP at the other interface and then couples 
to the outgoing radiative mode and exits. The transmit- 
tance is enhanced in this the photon can use two 
intermediate states (the first-order SP at both interfaces) 
to cross the metal film, but the enhancement is not as ef- 
ficient as when the hopping occurs back and forth several 
times building up a constructive interference. 

Absorption introduces another time into the problem: 
tabs, the typical time it takes for a photon to get ab- 
sorbed. Provided this time is smaller than both tres, and 
trad, the physical picture in not greatly altered. This is 
the case for the experimental parameters: absorption re- 
duces the Too maxima by a factor around 2-3, without 
altering either the position or the width of the peaks. 

The above discussion is for the special case of match- 
ing dielectrics (ei = £3) but the picture still applies in the 
non-symmetric case (ei £3), the situation in which ex- 
traordinary transmittance was initially discovered. Here, 
SP on the two different interfaces of the metal film have 
different frequencies and, instead of a SP diatomic ho- 
mopolar molecule forming in the presence of a continuum, 
we have instead a diatomic heteropolar molecule. From 
our model we find that there still are resonant peaks, 
but the transmittance maxima is reduced with respect 
to the symmetrical case. A detailed comparison between 
experimental and theoretical results in this case will be 
presented elsewhere. 

In summary, we have presented a detailed theoretical 
study of the extraordinary transmission of light through 
subwavelength hole arrays in metal films, obtaining good 
agreement with experimental data on free-standing films. 
An analytical minimal model, based on the numerical 
calculations, conclusively shows that the the enhanced 
transmission is due to tunneling through surface plas- 
mons. Different regimes in the enhanced transmission 
are found: for small film thickness the tunneling is reso- 
nant through "plasmon molecule" levels while, for large 
thicknesses, the photon hops from surface plasmon to sur- 
face plasmon, but exits the structure before the molecule 
level is developed, the tunneling being then sequential. 
Undoubtedly, this found ability of surface plasmons to 
transmit and focus light very efficiently is not restricted 
to this particular geometry and will be exploited for con- 
trolling light in new photonic devices. 
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FIG. 1. Experimental zero-order transmittancc of a square 
array of holes (lattice constant L = 750nm, average hole 
diameter of 280nm) in a free-standing Ag film (thickness 
h — 320nm). Inset: electron micrograph of the perforated 
metal film. 

FIG. 2. Calculated Too for an array of holes in a Ag film, 
defined by L = 750nm, h = 320nm and three different hole 
side lengths d. 

FIG. 3. Comparison between the calculated Too within the 
full (solid line) and minimal (dashed line) models, for the 
system considered in Fig.2, with d = 280nm. In panel (a) the 
dielectric constant is that of Ag, whereas in (b) absorption of 
Ag is neglected. 

FIG. 4. (a) \p\ versus A for the parameters defined in Fig. 
2, within the minimal model with no-absorption (solid line). 
Also e'''*'^''' for different values of h (in nm). (b) Too versus A 
for the metal thicknesses considered in panel (a) (successive 
curves being shifted by 100%). 
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